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The semiclassical propagation of spin coherent states is considered in complex phase space. For 
two time-independent systems we find the appropriate classical trajectories and show that their 
combined contributions are able to describe quantum interference with great accuracy. Not only 
the modulus but also the phase of the quantum propagator, both dynamical and geometric terms 
combined, are accurately reproduced. 

PACS numbers: 03.65.Sq, 03.65.Vf 



in 
o 
o 

(N 

■4— > 

o 
O 
o 

(N 

(N 
> 

(N 
(N 

in 
o 
in 
o 

S3 

» ■ 

G ■ 

s : 



I. INTRODUCTION 

The coherent state path integral for spin systems and 
its semiclassical approximation have appeared simulta- 
neously 0, and since then they have been inten- 
sively studied. Solari Q, and independently Kochetov 
|4|,y|, have shown the existence of an initially unexpected 
term that is sometimes called the Solari-Kochctov "ex- 
tra phase" , even though it is usually not a phase. Some 
important applications of spin path integral were in the 
study of spin tunnelling in the semiclassical limit 
even though initially some considered the method to be 
inaccurate Q . Stone et al E3 an d a l so Vieira and Sacra- 
mento have derived the spin coherent state semiclas- 
sical propagator in detail, paying particular attention to 
the Solari-Kochetov correction. This correction is related 
to the difference between the average value of the Hamil- 
tonian in coherent states and its Weyl symb ol 1 131 . and 
has a counterpart in the canonical case 0, Il3|. which 
has a flat phase space. The semiclassical description 
of tunne ling was reconsidered lately, using the instanton 
method HI EH- 

Recently, a semiclassical quantization condition for one 
dimensional spin systems has been derived |l6l | , including 
the first quantum corrections, in the spirit of the Bohr- 
Sommerfcld formalism. The same quantization condition 
was reobtained in [T3, where the authors also presented 
a semiclassical expression for the Husimi functions of sta- 
tionary states. Semiclassical theories for spin-orbit cou- 
pling in connection with trace formulas have appeared 
|18| and recently received renewed attention [l9j . 

It is interesting to note that the coherent state spin 
path integral is a natural setting for the investigation of 
geometric phases [2(J, |21j, a topic that has not only a 
fundamental importance in the mathematical structure 
of quantum theory |22| but also as an ingredient in the 
implementation of some quantum information protocols 
[23| . The semiclassical approximation to the spin propa- 
gator K(zi, Zf, T) = (zf\e~ lHT / h \zi) may sometimes lead 
to an expression to the geometric phase in terms only of 
classical quantities. 

In this work we present a concrete application of the 
semiclassical spin propagator. In general the calculation 
involves a classical trajectory, (z(t),z(t)), that starts at 
z — Zi and ends at z — z*j after a time T. Under these 



too stringent conditions, the only way to find a classical 
trajectory is by allowing the variable z(t) to be different 
from the complex conjugate of z(t) (which is denoted by 
z*(t)). We must therefore find a trajectory in C 2 that sat- 
isfies the boundary conditions z(0) = Zi and z(T) = z%. 
This is known as the root-search problem. Similar calcu- 
lations have already appeared for the canonical coherent 
states [24|, even in chaotic cases |25|. but so far no nu- 
merical example for a spin system has been presented. 

The article is divided as follows. In the next section 
the semiclassical theory of the spin coherent states prop- 
agator is briefly reviewed. The geometric phase and the 
tangent matrix associated with classical trajectories are 
also presented. In section III we study the simple sys- 
tem H = H 2 vJ 2 , which already has nontrivial proper- 
ties. In section IV the less symmetrical Hamiltonian 
H = h 2 v[J 2 + ir J 2 } is considered, and in section V we 
present the conclusions. 



II. THE SEMICLASSICAL PROPAGATOR 

Let \z) denote non- normalized spin coherent states, de- 
fined by 



\z) = exp{zJ+}| -j) = 



m=-j 



2j 



z j+m \m) 



(1) 



where the states \m) are the usual spin basis, and let 

K(z i ,z f ,T) = {zf\e- iftT / h \z i ), (2) 

be the quantum propagator, where H is the spin Hamil- 
tonian and T is the time. It has been shown that, in 
the semiclassical limit 7 — > 00, H = X/j — > 0, this can be 
approximated by 



Limit 7 

urn 



s K^f* elB/Rj d 2 s y r i \ 

(3) 

where the sum is over different classical trajectories. The 
exponent is the classical action S plus an extra term, the 
Solari-Kochetov (SK) correction: 



$ = S + 1 SK = S 



A{t)dt. 



(4) 
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It is well known that this semiclassical approximation is 
exact if the Hamiltonian belongs to the su(2) algebra 

The function A is given by 



A = 



d 



1 dH d 



1 dH 



dz4g(z,z) dz dz4g(z,z) dz ' 



(5) 



where H(z,z) = (z\H\z) / (z\z) plays the role of the clas- 
sical Hamiltonian and the metric factor is 



g{z,z) 



o 2 



■hi{z\z) 



2j 



dzdz 

The classical action is given by 



{l + zzf 



s = 



ihj- 



1 + zz 



H(z,z) 



dt + B, 



(6) 



(7) 



where B = —ihjln[(l + Zfz(T))(l + z(Q)zi)] is a boundary 
term, that takes into account the fact that in general z 
is not the complex conjugate of z. 

The integrals in Q and J7J are to be calculated along 
classical trajectories satisfying the Hamilton equations of 
motion 



1 dH 



h g(z, z) dz ' 
with boundary conditions 

z(0) = zt, z{T) = z). 
Notice that, in general, 

z(T)^z f , z(0)^z*. 



i 1 dH 

h g(z, z) dz ' 



(8) 



(9) 



(10) 



Since z(t) ^= z*(t) the Hamiltonian H(z,z), the action 
and the SK correction will in general all be complex num- 
bers. Therefore the term exp{ih~ 1 ^(zi, z/, T)} is usually 
not just a phase. 

If the propagator is written as a modulus times a phase, 
K(zi,Zf,T) = \K\e llp , it is well known that tp contains 

not just the dynamical term, —i Hdt, but also the ge- 
ometric (or Berry) phase ip g , which depends only on the 
geometry of the path traced by the state in the Hilbert 
space [2JJ, • This was shown by Berry in the context 
of Hamiltonians depending on slowly varying cyclic pa- 
rameters, and later generalized in many different ways, in 
particular to parameter independent, non-adiabatic and 
non-cyclic evolutions such as the one considered here [2l| . 
From equation @ we see that in the semiclassical limit 
the total phase will result from the interference of many 
classical trajectories, and that each one of them has an 
individual phase which is the sum of a dynamical term 
plus a geometric term given by 









/'I 

Jo . 



ihj- 



1 + zz 



A(t) 



dt 



(11) 



where ip p is a phase coming from the prefactor. 



If we set Zf = Zi in the simplest case H = vhJ z , we 
have z{t) = e+ ll/t Zi and z(t) = e w( *~ T) z*, and there is 
only one classical trajectory. The stability matrix ele- 
ment M zz (see below) in this case is simply e luT . For 
the very particular time T = 2'k/u it happens that 
z(t) = z*(t) and thus B is purely imaginary. If we use 
stereographic coordinates z = e 1 ^ tan(#/2) then it is easy 
to see that, since A and tp p both vanish, the geometric 
phase reduces to the well known result ip = 27r/'(l— cos 6). 
Generically, i.e. if H = vJ z but T is not a multiple of the 
period, or for more general Hamiltonians, the fact that 
z(t) 7^ z*(t) will introduce additional contributions com- 
ing from A, ip p and B, even for cyclic evolutions. Besides, 
more than one classical trajectory will be necessary, as 
we will see in the next sections. 

In the study of semiclassical spin tunnelling, the rele- 
vant classical trajectories are instantons with the remark- 
able property that z(0) = z* and z(T) = z } 14] . This 
happens because the initial and final points of the instan- 
ton minimize the average value of the Hamiltonian, and 
in that case the calculation is greatly simplified, leading 
to analytical results. However, as noted in |lq ] a more 
realistic description of the system requires the addition 
of higher order terms to the Hamiltonian that destroy 
this simple property and make it necessary to consider 
the more generic trajectories we have discussed (which 
have been called "boundary jump instantons" in [la])- 

Note that the action (7J , with the necessary boundary 
term, leads to the following Hamilton- Jacobi relations 



i dS _ 2jz(t) ^dS_ 
hdzj l + zfz(t) hdzi 

dS H 



2j*(0) 
l + *(0>i' 



(12) 
(13) 



Note also that if the Hamiltonian is 0(j), then S is O(j), 
but the SK correction is 0(1), and therefore can be con- 
sidered small in the semiclassical limit. 

The prefactor is related to the tangent matrix, or sta- 
bility matrix, as follows. Small variations in the bound- 
ary points Szi and 5zf induce variations Sz(T) and 5z(0). 
Taking derivatives of the Hamilton- Jacobi relations 112|) 
it is possible to write 



(l + z f z(T)r 2 6z(T)\_ 
(l + z t z(0))- 2 8z(0) 



A A - 

Azz Azz 



Szi 
Sz f 



(14) 

where all matrix elements can be written in terms of sec- 
ond derivatives of the action (see [lj|] for an analogous 
calculation with canonical coherent states). On the other 
hand, the tangent matrix of a given trajectory (z(t), z(t)) 
is defined as the linear application that takes a small ini- 
tial displacement to a final displacement, 



Sz(T) 
Sz f 



M zz 
M- Z2 



M zz 



Szi 
Sz{0) 



(15) 



Manipulating equation i|14[l one can show that the ele- 
ment M zS is related to the prefactor in © according to 



3 



d 2 S 



2.7 



hdz.dzf (I + z t z(0)) 2 M sg ' 



(16) 



Since the tangent matrix may be numerically integrated 
together with the coordinates, this expression is very con- 
venient in practice. The phase tp p may then be followed 
dynamically, imposing that K = 1 (ip p = 0) for T = 0. 

In practice, one must find all values of z(0) for which 
z(T) = z%. Notice that these initial values are implicit 
functions of the time T, and as such they trace out curves, 
or "branches" , in the z(0)-plane. There will in general 
exist more than one branch for each value of T, and one 
must add all contributions coherently. Note, however, 
that the contribution of a branch may display an erro- 
neous increase and need to be removed for a certain re- 
gion of the parameter T, because of a phenomenon that 
is common in as ymp totic expansions known as Stokes' 
phenomenon |26l 1271 1281 |29| . A discussion of this prob- 
lem in the context of the one-dimensional semiclassical 
propagator may be found in |3p| . Stokes' phenomenon 
will not be relevant to the present work. 

From equation (|16[) it is clear that the semiclassical 
approximation fails in the vicinity of points for which 
MzzjT) = 0, which are called phase-space caustics 
|25l l30l l3l| . Semiclassical uniform approximations that 
are valid near caustics have been obtained for the canon- 
ical coherent state propagator using a conjugate of the 
Bargmann representation in |32j |. and a similar calcula- 
tion is in principle possible for the spin propagator, but 
here we shall not be concerned with the effect of caustics. 



III. FIRST EXAMPLE, H = h 2 vJ 2 z 



We consider a simple Hamiltonian, 
H = h 2 vj 2 



(17) 



(the constant v has the appropriate units) and we shall 
be interested only in the diagonal propagator 



K{z u T) = ( Zl \e- lHT l h \z t ), 



(18) 



whose squared modulus, after proper normalization, cor- 
responds to the return probability as a function of time. 
This is given by 



j 

K( Zi ,T)= ]T \(m\zi)\ 
m=—j 



2 -imrhvT 



(19) 



which for integer j is periodic with period T r = 2n/hv. 
In the semiclassical limit the term that contributes the 
most to this sum (see JTJ) is m = j(l z i| 2 — l)/(|zi| 2 + 1). 
If we linearize the exponent in the vicinity of this term 
we have 



K{z u t) 



(m \z t )\ 2 e- im ° T e- 2m ° nT , (20) 

n~mo 



where we have introduced a scaled time 
t = hvT. 



(21) 



Notice that expression Il2()f has a different time scale, 
t c = n/mo. The quantities T r and T c = T r /hu are usually 
called revival time and classical time [33j. 

Let us turn to the semiclassical approximation, which 
has been analyzed in some previous works 0, llOl lll| . 
The classical Hamiltonian is 



H = n 2 uj(j-l] 



i 



zz ■ 



1 



h 2 uj 



(22) 



The classical equations of motion, z = —ihvujz and z = 
ihvuiz (here a dot denotes derivative with respect to t), 
together with the boundary conditions z(0) = Zj and 
z(T) = z* , have the simple solutions 



z(t) 
Notice that 



z(0), z(t) = e 



2/.t 
3 



lhuuJt z(0). 



zz+1 



(23) 



(24) 



is a constant of the motion, where fj, = j(j — 1/2). Cal- 
culating it at the initial and final points we have the 
consistency condition 



2/j, 
3 



M 2 



i 



+ i 



(25) 



which in general has an infinite number of solutions. We 
come back to that later. Notice also that for z(0) = z% 
and z(0) = z* the motion is periodic with period 2n /hvu, 
which in the semiclassical limit becomes T c . 
The action can be easily found to be 



- = -2yIn(l- 



■e- iUT \zA 2 ) 



4 M 



and the Solari-Kochetov term is also available, 

,2 " 



A = W-v 



2j 



8/u 



(26) 



(27) 



To find the prefactor, we consider small variations Szi 
and 5z(0). The final value of z will be changed according 
to 

z{t) + 5z(t) = e l ^ +s ^ T (z(0) + 5z(0)) 

* z{t) + M sg 8z(0), (28) 



which leads to 



4/i irziz(0) 
J(l + z l z(0)) 5 



(29) 



There are two different phase space caustics, that can 
be found by equating = 0. They are located, as 
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FIG. 1: The normalized propagator as a function of scaled 
time r = KvT (in units of the classical period r c ) for Zi = 0.5 
and j = 50. In (a) we see the separate contributions of 6 
branches (each peak is due to a different branch) . When they 
are added coherently we get (b), which is indistinguishable 
from the exact result, including the fast oscillations. 



functions of the scaled time, along the curves z±(0) = 
— z^ 1 (l + ir =F \J2it — t 2 ) . We shall not be concerned 
here with the effect of caustics. 

The problem now is to find the values of z(0) for which 
z(T) = z* . These points form curves in the plane f(0), 
that we call branches and denote 2 T (0), since they are 
parametrized by the time. We note again that even for a 
fixed time there may exist many branches. The method 
used for finding them was the following: for a fixed value 
of time, To, a regular grid is placed on the z(0) plane, and 
each point is taken as an initial condition. Those points 
for which z(tq) is close to z* are then used as initial 
guesses for a root-finding procedure. For r = To ± St the 
previously obtained solutions serve as initial guesses, and 
no grid is used. This way the branches may be obtained 
rapidly and exhaustively. Of course the point z(0) = z* 
is a solution for r = 0, and thus it must be contained in 
one of the branches. 

As an example, let us consider Zi = 0.5 and a very 
semiclassical regime, j = 50 (which gives h = 0.02). For 
the time interval of five classical periods we have found 



FIG. 2: The different branches in the z(0) — x + iy plane 
contributing to the semiclassical return probability, for the 
same parameters as in the previous figure (quantities are di- 
mensionless) . These curves as parametrized by the time, and 
each one of them passes through the point (0.5, 0) at a differ- 
ent instant. Only the vicinity of this point is shown, for the 
sake of visibility. 



6 different branches. The squared modulus of each of 
their normalized contributions is shown in Fig. la, with 
different line styles. We see that each peak in the return 
probability is due to a different branch. The appearance 
of the classical time scale t c w 0.1, which comes from 
mo = —30, is very clear. For r « the branch depicted 
in solid line is in the vicinity of the point z*, and as r 
approaches the integer multiples of r c each one of the 
other branches approaches this point, generating a peak 
in the return probability. The behavior of the different 
branches can be seen in Fig. 2. After the third peak 
the contributions start to overlap, and therefore we must 
add them coherently to get interference effects. When 
this is done, the result is as shown in Fig. lb. The exact 
calculation is easy to perform, and it is indistinguishable 
from the semiclassical one at this scale. Even the fast 
oscillations for r > 2.5r c are reproduced. 

When the value of Zi corresponds to a point along the 
equator of the Bloch sphere, i.e. when \zi\ — 1, the calcu- 
lation may be simplified. In that case there exists a very 
particular classical trajectory: the one for which z(t) = Zi 
and z(t) = z* . In this case it easy to see that u> = and 
there is no movement. Let us call this the static trajec- 
tory. Each different value of z T (0) determines a certain 
ui, and the final expression for the normalized propagator 
is 



K sc {Zi,T) 



exp 



IT 


(- 2 




\irj ( 


'lu 2 




2j 



2„-2 



V) 



;2 



(30) 



where we have used (|24f) and the fact that zz is a constant 
in time. It is easy to see that each term in this sum (as 
well as equation l|25|) ) is even in u>, and therefore that 
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FIG. 3: The normalized propagator as a function of scaled 
time (in units of r r ) for z% = 1 and j — 10. In (a) we see 
the separate contributions of 11 branches (all of them except 
one must be counted twice, see text). When they are added 
coherently we get the solid line in (b) . The dashed line is the 
exact result. 

trajectories must come in pairs. The net effect is that 
we may search only for a/s with a positive real part, and 
once a certain trajectory is found its contribution must 
be doubled, except for the static one, which gives 

\K st ( Zl ,r)\ 2 = fl + r^j-l)^ \ (31) 
For short times we may approximate 

|X st (^)| 2 «l-y(j-^) 2 , (32) 

and comparing this with the exact short-time return 
probability 

\K(z u r)\ 2 ^\-{{H 2 )-{H) 2 )T 2 /h 2 

= l-\j{j-\), (33) 

(where (■) = • |zj)) we see that in the semiclassical 
limit the short-time regime is well described by the static 
trajectory alone. 



FIG. 4: The semiclassical (solid) and the exact (dashed) cal- 
culations of the real part of the normalized propagator, for 
Zi = 1 and j — 10. The imaginary part is also reproduced 
with great accuracy. Therefore the total phase, dynamical 
plus geometrical terms, can be recovered from the semiclassi- 
cal expressions. 

As a concrete example, we have considered zi = 1 and 
j = 10, so that higher order quantum effects may become 
more visible. In this case 11 branches were found that 
contributed to the final result. Their separate contribu- 
tions to the normalized propagator can be seen in Fig. 
3a as functions of time (in units of the revival time r r ). 
Because of the even parity in u>, each one of them was 
considered twice, except for the static one. When we add 
them coherently, we get the solid line in Fig. 3b, which 
is to be compared with the exact result, the dashed line. 
We see that the later is reproduced with an extraordinary 
accuracy, including the quantum revival. 

We do not show the branches in the z(0) plane, but 
their behavior is quite simple. The static branch lyes in 
the point z(0) = 1 for all times. For short times all other 
branches are at the vicinity of — 1 , but their contribution 
to the propagator is negligible. As time passes, they move 
away from this point and towards 1. The most impor- 
tant thing to note is that all branches have appreciable 
contributions at the period t = 2n, so the phase of the 
propagator is determined by their coherent superposition 
and has no simple expression. 

Concerning the phase in (|30() . notice that it can be 
written as %t{w 2 — l)/4 because fi/j 2 — 1 is actually equal 
to —l/2j. Since r is of order h the second term vanishes 
in the semiclassical limit, but the first one remains be- 
cause u> 2 ~ j 2 . However, we have observed numerically 
that the factor e~ lr / 4 , which is common to all branches 
(and thus has no relevance to the modulus), destroys the 
correspondence with the exact result. Once this term is 
removed, which is equivalent to the semiclassically ac- 
ceptable approximation 

M = j(i-l/2)«j 2 , (34) 

the agreement is excellent, as we can see from Fig. 4, 
where the real part of the normalized propagator is shown 
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FIG. 5: The normalized propagator as a function of scaled 
time (in arbitrary units) for Zi = 1 + i and j = 10 in the 
second example. In (a) we see the separate contributions of 6 
branches. When they are added coherently we get the solid 
line in (b). The dashed line is the exact result. 



(the imaginary part is also very accurate) . This approxi- 
mation must not be done in the prefactor, where we have 
used n — j(j - 1/2). 

We therefore conclude that expression ® is able to 
reproduce both modulus and total phase, dynamical and 
geometric terms combined, of the quantum propagator 
with great accuracy in the semiclassical limit. It involves 
a coherent sum over classical trajectories, and thus it 
is not possible to extract the geometric phase from it 
in a closed form. As noted in the previous section, the 
geometric phase will be given by (|llfl only in the simplest 
cases, when a single classical trajectory is involved. Even 
then it will in general be different from the simple formula 
ijdt(zz — zz)/(l + zz). 

Stone et al have analyzed the case H — J 2 in detail 
|lf)j | . They have shown the formal equivalence of the semi- 
classical and exact (non-normalized) propagators up to 
order j°, and also that for \zi\ = 1 the result for a mas- 
sive particle constrained to move on a ring is recovered 
as j — > oo. However, in |lf| no numerical calculation has 
been done, and the fact that replacing [i by j 2 in the 
phase improves the result was not noticed. 




FIG. 6: The semiclassical (solid) and the exact (dashed) cal- 
culations of the real part of the normalized propagator, in the 
same situation as the previous figure. Once again the phase 
is accurately reproduced. 

IV. SECOND EXAMPLE, H = h 2 u[J 2 + irj 2 ] 

As a second example, we choose the less symmetrical 
Hamiltonian H = h 2 v[J 2 + kJ 2 ], which has been con- 
sidered in 0, ^3 • We have chosen the irrational value 
k = it, in order to avoid too much simplicity. The exact 
calculation of the return probability can be done by di- 
rect diagonalization of the Hamiltonian in the J z basis. 
The initial state was taken as a complex number out of 
the equator, z; = 1 + i. For very large values of j the 
return probability is very similar to what we see in Fig. 
lb, but for more moderate values it is much less regular. 

In Fig. 5a we see the separate contribution of 6 
branches, each of which has only one peak. This time 
the classical trajectories, the tangent matrix, the action 
and the SK correction all had to be determined through 
a numerical integration, using a Runge-Kutta method. 
The coherent superposition of the individual contribu- 
tions is the solid line in Fig. 5b, and the dashed line is the 
exact result. Even though the dynamics does not show 
any periodicity or regularity in the considered time in- 
terval, the semiclassical approximation works very well. 
Notice that we have used the same scaled time as in the 
previous example, t = hvT. 

The phase of the propagator is again recovered with 
great precision, as we see in Fig. 6, where the real part of 
both the semiclassical and the exact calculations is shown 
(again the imaginary part is reproduced with the same 
accuracy). This time no approximation of the kind l|34|) 
was made to the phase, which was obtained numerically. 



V. CONCLUSIONS 

A numerical evaluation of the semiclassical spin coher- 
ent state propagator has been presented for the first time. 
We have considered two different systems as examples, 
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the quantum rotor H = h 2 vJ 2 and the less symmetrical 
Hamiltonian H = Y?v\J 2 + irJ 2 ]. In the first case the 
exact calculation of the return probability is very simple, 
and is reproduced by the semiclassical one with excellent 
precision for j = 50. For the smaller value j = 10 the 
accuracy of the approximation is still remarkable. 

We have made the numerical observation that a factor 
—ihvT/4 must be removed from the semiclassical phase 
in order to reproduce the exact quantum phase for the 
rotor (the modulus is not affected) . We have no rigorous 
explanation for this fact, except that this term vanishes in 
the semiclassical limit and that its removal corresponds 
to the replacement j(j — 1/2) — > j 2 . Since phases are 
defined modulo 2ir, we believe this limit must be con- 
sidered with great care. Since no such adjustment had 
to be made for the second Hamiltonian studied, we con- 
clude that geometric phases are adequately reproduced 
by semiclassical calculations. However, they generally 
do not have a simple closed formula, but result from the 
interference of many classical paths. 

The most important difficulty in the implementation 
of the semiclassical spin propagator in practice is finding 
the relevant classical trajectories, i.e. those values of 1(0) 
for which the boundary condition at time T is satisfied, 
z(T) = z*f. We have done this by first propagating for 
a time T a grid in the z(0) plane, thus finding initial 
rough estimates for the relevant points. These were then 
used to feed a root-finding procedure. The whole process 
is rather artisanal and hard to automatize, making it 
almost impossible to tackle problems in which a large 



numbers of trajectories is necessary. The same kind of 
difficulty is found in the case of canonical coherent states. 

Another problem, which we have not considered, is 
the existence of caustics, initial conditions for which the 
prefactor diverges at time T. We have seen that for 
H = h 2 vJ 2 there are two such points, but they did not 
have any influence in the cases we have analyzed. For 
systems with higher nonlinearities it is reasonable to ex- 
pect a larger number of caustics, which could also hinder 
the practical application of the semiclassical approxima- 
tion. In (3^ the authors have developed an approach to 
treat the vicinity of caustics in the canonical coherent 
states propagator, that makes use of a certain dual of 
the Bargmann representation. This could in principle be 
generalized to the spin case. 

The main virtue of the semiclassical propagator is its 
weak dependence on the dimension of the representation, 
2j + 1. If this is too large a direct diagonalization of 
the Hamiltonian may become impracticable, while the 
difficulty with the semiclassical method is basically the 
same. Therefore an excellent approximation to the exact 
result, at least for short times, would be easily available. 
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